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Density-functional calculations are carried out to investigate incommensurate magnetic structures and ferro- 
electric polarization in newly discovered multiferrioc material Mnl2 . We calculate magnetic exchange coupling 
parameters and construct an effective magnetic model whose magnetic orders are in good agreement with those 
measured experimentally. The measured direction of polarization is also consistent with the result derived from 
the symmetry analysis of the magnetically ordered state. In particular, we find that the inter-plane magnetic 
exchange coupling is pivotal to the emergence of the spiral magnetic structure that induces ferroelectricity. 
Our calculation suggests that spin-orbit coupling on I ions is primarily responsible for the appearance of the 
polarization. 

PACS numbers: 75.85.+t, 75.10.Hk, 71.70.Ej, 71.15.Mb 



I. INTRODUCTION 

Multiferroics, which exhibit magnetic and dielectric or- 
der in the same phase, recently have attracted increasing 
attention^'^. The recent experimental research on multifer- 
roics has shown that ferroelectricity (FE) and magnetism cou- 
ple so strongly that the electric degree of freedom can be ma- 
nipulated by an external magnetic field and vice versa^"^^. 
These properties offer unprecedented applications in modem 
energy-effective electronic data storage technology 

Theoretically, phenomenological models and symmetry 
analysis have clarified the circumstances where a spiral spin 
structure can induce an electric polarization^ Harris^^ 
gives a simple method to describe the magnetic ordering and 
their relationship to ferroelectricity based on lattice, space and 
time reversal symmetries: the symmetry of the magnetoelec- 
tric interaction can determine the direction of the magnetic 
induced spontaneous polarization. 

Several microscopic mechanisms have been proposed to 
explain the magnetoelectric coupling in multiferroics. One 
is the well-known Katsura-Nagaosa-Balatsky (KNB) model 
which is based on the idea that spin currents are induced 
between the spiral spins and can therefore be considered 
as electric moments. The second is that magnetically- 
induced ionic displacements induce the polarization through 
Dzyaloshinskii-Moriya (DM) interactions The electric 
cancelation model^^ gives a simple but general mechanism to 
understand the interplay between ferroelectricity and nonlin- 
ear magnetism in multiferroics. As a powerful tool to inves- 
tigate the electronic structure of materials, density functional 
theory (DFT) has played an important role in the understand- 
ing of the collinear-spin type^^'^^ and the spiral magnetic ma- 
terials LiCu202 and LiCuV04^^ 

Mnl2 has been investigated primarily from the interest in 
magnetic and optical properties^"^"^^. However, it was discov- 
ered recently by Kurumaji et al?^ that Mnl2 is also a mul- 
tiferroic material. Mnl2 crystallizes in the Cdl2 type struc- 
ture with the space group P3ml (No. 164). The unit cell con- 
tains one molecule with the manganese ion located at (0, 0, 0) 



and the iodide ions at |, ix), where u = 0.245 ± 0.002, 
a = 4.146 A and c = 6.829 A^"^. Magnetic properties are 
dominated by Mn^^ ion with S = \. Sato et alP ob- 
served three successive phase transitions at 3.95 K (Tati), 3.8 
K {Tn2) and 3.45 K (T/va). As temperature is decreased, 
the Bragg reflection at q^^ (0.1025,0.1025,0.5) appears 
at Tni- When the temperature is further decreased, the re- 
flection position begins to move slightly out of the {hhl) 
plane towards the {h{)l) plane. Finally at T/v3 it jumps to 
q^^ (0.181,0,0.439), in which we notice that hi is not 
equal to /c2, where ki is the i\h component of the magnetic 
modulation vector. Blow T/v3, the proper screw magnetic 
structure is realized, which induces FE polarization about 84 
liC /m? along [110] direction at 2 K.^^ Moreover, an in-plane 
H can induce the rearrangement of the six multiferroic do- 
mains and every 60° rotation of the in-plane H leads to 120° 
flop of the P direction as a result of the flop of magnetic order. 

Important questions concerning Mnl2 are why it has the 
helix spin magnetic ground state and how the spiral spin in- 
duces ferroelectric polarization. It is also of great interest in 
the appearance of successive phase transitions as temperature 
decreases. In this paper we perform a comprehensive theo- 
retical investigation of these intriguing properties. We first 
calculate the magnetic exchange coupling parameters in Mnl2 
and then discuss the phase transitions within mean-field the- 
ory based on a Heisenberg-type magnetic exchange Hamil- 
tonian, in which six exchange interactions are taken into ac- 
count. Then we calculate the polarization of Mnl2 and per- 
form symmetry analysis to show the polarization is consistent 
with magnetic order. We find that the inter-plane coupling is 
fairly strong because of linear exchange path and is extremely 
important in inducing the spiral magnetic order ground state. 
Finally, we show that the spin-orbital coupling (SOC) on I 
ions makes primary contribution to FE polarization, based on 
an analysis of charge density difference between cases with 
SOC and without SOC. 

The paper is organized as follows. First, in Sec. II, we per- 
form DFT calculation to obtain the six exchange parameters 
from seven spin ordered arrangements and determine the mag- 
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netic vector of Mnl2 by using these exchange parameters. 
Then, in Sec. Ill, we determine the direction of the polariza- 
tion in Mnl2 through symmetry analysis and calculate the po- 
larization of different magnetic vectors using DFT. Finally, in 
Sec. IV, we give a summary and provide the main conclusions 
of our paper. 

II. WAVE-VECTOR SELECTION IN Mnis 

A. Calculation of the exchange parameter 

Our DFT calculations employ the projector augmented 
(PAW) wave method encoded in Vienna ab initio simula- 
tion (VASP) package^^"^^, and the generalized-gradient ap- 
proximation (GGA) for the exchange correlation functional^ ^ 
Throughout this work, the cutoff energy of 400 eV is used 
when expanding the wave functions into plane- wave basis. A 
set of 2 X 4 X 2 F centered k-points is used for the 4 x 2 x 2 
supercell in the calculation of exchange parameters, which 
is sufficient to obtain converged results for all quanties un- 
der consideration. The GGA plus on-site repulsion U method 
(GGA+U) in the formulation of Dudarev et al?^ is employed 
to properly describe the electron correlation effect associated 
with the Mn 3d states by an effective parameter Ueff.^^ Var- 
ious Ueff values for Mn are taken to obtain the exchange 
coupling parameters. The self-consistent-field convergence is 
achieved when the total electronic energy difference between 
last two cycles is less than 10~^ eV. In all our calculations, we 
use the experimental crystal structure directly taken from the 
experiment^"^. 

In order to obtain the exchange parameters from DFT calcu- 
lations, we separate the total energy into nonmagnetic(i^non) 
and magnetic contributions 




FIG. 1: (color online). Side view of the Mnl2 and the red and green 
atoms are Mn and I, respectively. Additionally, the exchange parame- 
ters Ji, Jab, J2, Jc, Jnc and Jnnc between the cations connected 
by arrows are defined. 

The six spin exchange parameters can be evaluated by ex- 
amining the eight ordered spin states of Mnl2 (i.e. the FM, 
AFMl, AFM2, AFM3,AFM4, AFM5, AFM6 and AFM7 
states) defined in Fig. 2 in terms of a 4 x 2 x 2 supercell. Table 
I summarizes the relative energies of these states per 4x2x2 
supercell, i.e. per 16 f.u., determined from our GGA+U and 
spin-orbit coupling (SOC) calculation. Then, by applying the 
energy expressions obtained for spin dimers with N unpaired 
spins per spin site (in present case, N=5)^^, the energies of the 
eight magnetic states per f.u. are written as 



H — H^Q^ + / ^ JijSi • Sj^ (1) 

i<j 

where Si and Sj are the spin operators on the spin sites i and 
j, respectively and the Jij is the exchange parameter between 
the sites i and j. Jij < corresponds to the ferromagnetic 
(FM) coupling between the two sites while Jij > corre- 
sponds to antiferromagnetic (AFM) interaction. 

Fig. 1 illustrates the magnetic pair exchange interaction 
used in our modeling. Ji, Jab and J2 are the intraplane in- 
teractions between the cations. Jc, Jnc and Jnnc are the 
interplane couplings between Mn^+. As it can be seen from 
the Fig. 1, the distance between the two Mn cations in the 
coupling Jnnc is much longer than that of the coupling Jn 
and Jnc, so we may expect that the coupling Jnnc is much 
weaker due to the decreasing of the overlap of the wave- 
functions atoms with the increasing distance between the two 
ions. However, according to Wollan et al?^ the interaction 
Jn nc is fairly strong since it has an almost linear exchange 
path (Mn^+-Br~-Br~-Mn^+) in MnBr2, whose structure is 
isomorphous with that of Mnl2 . Our calculation shows that 
the magnitude of Jnnc is almost the same as Jn and Jnc, 
which confirms their conclusion. 



Efm = ^(3Ji + 3J2 + 3Jab + Jc + SJatc + 3Jnnc), 

El = — Jl + 3J2 — Ja6 + Jc — ^Jatc — JatATc), 

E2 = — (3Ji + 3J2 + 3Jab — Jc — QJnc — 3Jnnc), 

E3 = -^{Jl — J2 — Jab Jc '^JnC — JnNc)i 

Ea = Jl + 3J2 - Ja6 - Jc + 2JArc + JatATc), 

E^ = ~^{Ji — J2 - Jab — Jc — '^Jnc + Jnnc), 
N'^ 1 1 11 

E7 = — - Jl + J2 - - Ja6 + Jc - ^"^^C' + 2"^^^^^' 

Eg = ^(Ji + J2 + Jc + 2J^c). (2) 



By mapping these relative energies onto the corresponding 
relative energies from the magnetic energies (spin-exchange 
energies, relative energies with respect to FM state), we ob- 
tain seven equations. But there are only six spin-exchange 
parameters, Ji, Jab, J2 ,Jc, Jnc and Jnnc- For this overde- 
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FIG. 2: (color online). Schematic plots of eight different magnetic orderings in Mnl2 used for GGA+U calculation to extract the five spin- 
exchange parameter Ji, J2, Jab, Jc, Jnc- The red and green circles represent the up and down magnetic moment on Mn sites, respectively. 



TABLE I: Relative energies (in meV) of the six ordered spin states of Mnl2 obtained from the GGA+U and with SOC calculation. In calculation 
with SOC, the spin is parallel with the c axis. The value in brackets is the Uef / in unit of eV. 



state 


GGA+U(0) 


GGA-hU(0)-hSOC 


GGA-hU(l) 


GGA-hU(2) 


GGA-hU(3) 


GGA-hU(4) 


GGA-hU(4)-hS0C 


FM 























AFMl 


-258.7 


-250.2 


-149.5 


-81.1 


-35.8 


-4.8 


-2.8 


AFM2 


-83.1 


-82.9 


-58.4 


-41.5 


-29.6 


-21.0 


-21.0 


AFM3 


-203.8 


-199.3 


-127.9 


-79.0 


-45.8 


-22.7 


-21.6 


AFM4 


-241.8 


-233.2 


-137.7 


-73.1 


-30.6 


-1.7 


0.3 


AFM5 


-210.7 


-206.4 


-133.2 


-83.4 


-49.8 


-26.3 


-25.3 


AFM6 


-166.6 


-162.3 


-101.3 


-59.7 


-31.9 


-12.6 


-11.5 


AFMl 


-253.6 


-245.6 


-150.4 


-85.3 


-41.9 


-12.0 


-10.2 



termined system of equations, we obtain these parameters by 
using a least-squares technique^^'^^, as summarized in Table 
II. 

It is noted from the Table II that the intraplane exchange 
couplings Ji , J2 and J ah are antiferromagnetic for Uef / < 3 
eV. However, Ji becomes negative (i.e. to be ferromagnetic 
coupling) with Ueff — 4 eV, which is consistent with the es- 
timated coupling in a similar compound MnBr2^^ in neutron 
diffraction experiment. As a consequence, the intraplane spin- 
exchange interactions are geometrically frustrated. We no- 
tice that the exchange coupling is rather weak with respect to 
similar compound such as CUCI2 (exchange parameter about 
10 meV)^^. The weak exchange coupling can also be seen 
from the low magnetic phase transition temperature (3.45K - 
3.95K). The intraplane exchange coupling Ji arises via direct 
exchange and superexchange interactions, mediated by two 
Mn-I-Mn bond with the same bondangle 90.44°. For the sim- 
ilar case of Cu-O-Cu bonds, it has been shown that when the 
bond angle is close to 90° the resulting exchange energy is 
rather small."^^ For the case of Uef / = 4 eV, one may notice 
that the coupling Ji becomes ferromagnetic and weaker than 
J2, which has the twice of the distance of Ji, so it seems really 



strange. This is because both direct exchange (ferromagnetic) 
and superexchange coupling contribute to Ji , which have the 
opposite sign, thus there is a cancelation. As the increasing 
of U, the superexchange coupling becomes weaker while the 
direct coupling is almost unchanged, which finally leads to 
a weak ferromagnetic coupling. Furthermore, the coupling 
is weaker compared to the antiferromagnetic coupling J2 be- 
cause of the cancelation. 

The interplane coupling parameters Jc, Jnc and Jnnc 
are AFM coupling for all Ueff, which is consistent with 
the experiment carried by Cable et al?^'^^ . The fact that 
this two-anion indirect exchange coupling Jn nc (mentioned 
above, has linear exchange path) appears to be antiferromag- 
netic might be expected by analogy with the single-anion su- 
perexchange mechanism. This coupling, which has linear ex- 
change path (Mn^+-I~-I~-Mn^+), is stronger with respect to 
the other two interplane coupling for Uef / < 3 eV. The mag- 
nitude of the three interplane interactions are almost the same 
for Uef / = 4 eV. 

From the theory of superexchange it follows that the cor- 
responding coupling constant is proportional to 1/U. So if a 
coupling is mediated mainly by the superexchange interact- 
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TABLE II: Values of the spin-exchange parameters J (in meV) in Mnl2 obtained from the GGA+U and GGA+U+SOC calculation. In calcula- 
tion with SOC, the spin is parallel with the c axis of the Mnl2. The value in brackets is the Uef / in unit of eV. 



J 


GGA-hU(O) 


GGA-hU(0)-hSOC 


GGA-hU(l) 


GGA-hU(2) 


GGA-hU(3) 


GGA-hU(4) 


GGA-hU(4)-hS0C 


Jl 


0.49 


0.46 


0.27 


0.13 


0.03 


-0.03 


-0.03 


J2 


0.14 


0.13 


0.10 


0.07 


0.05 


0.04 


0.04 


Jab 


0.03 


0.03 


0.02 


0.01 


0.01 


0.01 


0.01 


Jc 


0.04 


0.04 


0.03 


0.02 


0.02 


0.01 


0.01 


Jnc 


0.03 


0.04 


0.02 


0.02 


0.01 


0.01 


0.01 


Jnnc 


0.06 


0.05 


0.04 


0.03 


0.02 


0.01 


0.01 



ing, one will expect a strong influence of the Hubbard param- 
eter on the strength of this coupling. That is the reason why 
most of the exchange coupling strength significantly decrease 
with increasing Ueff- The variation of Jnnc seems to show a 
1/U dependence, so it is mediated by superexchange interac- 
tion, which is an apparent result by the exchange path analysis 
above. While Uef / does not have much influence on 7^6, thus 
it is a direct exchange coupling. The exchange parameter does 
not have much difference whether SOC is considered or not. 

Thus in Sec. II B, we use the coupling obtained from cal- 
culation without SOC. 



B. The classical groundstate of the Mnl2 

To make the problem not complicated, we will describe the 
magnetic ordering by a version of mean-field theory in which 
one writes the magnetic free energy Fm as 

^=\Y. X-l{r^,r,)Sa{r.)Sp{rj) + 0{S% (3) 

Yi^a]Yj ,/3 

Where Sai^i) is the thermally averaged a component of 
the spin at position r^. We introduce Fourier transformations 
of Sa (fi ) and (^i -,^3) ^ we omit high order terms thus have 



SM.n) = 1^5„(R + Ti)e^''-(«+-*), (5) 



R 



1 

N 

R 

where N is the number of the unit cells in the system, is 
the location of the ith site within the unit cell, and R is the lat- 
tice vector. As our main interest lies in explaining wave- vector 
selection, we have completely ignored anisotropy, whose ma- 
jor effect is to select the spin orientations. That is an isotropic 
model, so that we have 



(7) 



is the entropy associated with a spin S. In our case, we 
only consider the exchange coupling defined in Fig. 1. In 
our coordinate system, the lattice vector is di = ai, a2 = 
— \ai + ^aj and as = ck (see Fig. 6, xl, yl), where a and 
c is the lattice constant of Mnl2 . Thus the reciprocal lattice 
is hi = ^(?+ ^J), h2 = ^ • ^Jand 63 = ^k. We 
set q = kihi ^ k2h2 ^ k^bs in the Fourier transformation, we 
have the following x~^{q) with some algebra in Mnl2, 

X-\q} = K^ckT 

+ 2Ji[cos(/ci) + cos(/ci + /C2) + cos(/c2)] 

+ 2Ja5[cos(2A:i + A:2) + cos(A:i +2^:2) 

+ cos(-/ci + /C2)] 

+ 2J2[cos(2fei) + cos(2A;i + 2^2) + 008(2^2)] 

+ 4JArc[cos(/ci) + cos(/ci + /C2'* ^ '"'^'^''^ 



COs(/C2)] COs(/C3) 



where 6^(3 is unity if a = (3 and is zero otherwise and is a 
spin-dependent constant of the unity, so that —ck 5'a(r)^ 



2Jc cos{ks) + 2 JArArc[cos(/ci + 2/c2 + ^3) 
+ cos(2A:i + A:2 - ks) + cos(A:i - A:2 + ks)]. (8) 

We set x~^{q) ^ + ckT + J{q) and substitute 
the exchange parameters calculated with Ueff = 4 eV 
(GGA+U4+S0C has little effect on our results) in Table II 
into Eq. (8). We also use data from other Ueff and find 
that the magnetic vectors have much difference with those 
in experiment, so here we just list the result of Uef f = 4 
eV and we find that the magnetic moment of Mn calculated 
with Ueff = 4 eV is 4.58 /i^, which is very close to the 
experimental moment^"^. To explain different q at transition 
point in experiment, we need to find the value of q to min- 
imize X~^(q) according to q selection^^. At T/vi and Tns, 
the picture of X~^(^i7 ^2) is in Fig. 3 and Fig.4, respectively. 
At Tni (we set /c3 = 0.5 ), from Fig. 3, the minimum 
point of x"^ is at {ki, k2) = (0.1226,0.1226), which is in 
good agreement with experiment value (0.1025, 0.1025). At 
Tn3 (set ks = 0.439), from Fig.4, minimizing x~^ we get 
(/ci, /C2) = (0.155,0.089), which is also consistent with the 
experimental value (0.181, 0). In Fig. 5, we fix k2=0 and min- 
imize x~^' then we obtain (/ci, ks) = (0.206, 0.444), which 
is close to (0.181, 0.439) below Tns in experiment. 

It is of interest to notice that ki and /c2 become asymmetric 
when temperature is below Tns while ki equals to k2 at Tni- 
When the value ks is special (0 or 0.5), J{q) is invariant for 
the exchange ki and /c2, thus we will expect ki and k2 have 
the same value. However, if ks is not special (i.e. neither 
nor 0.5), from Eq.8 we find that only the coupling Jnnc have 



FIG. 3: (color online). The diagram of (/ci, A;2, = 0.5) 
at Tni- It is noted from above picture the minimum point is at 
{ki, k2) — (0.1226, 0.1226), which is close to the experiment value 
(0.1025,0.1025). 



FIG. 6: (color online). The coordinate system in the paper, a denotes 
the lattice basis while b denotes the reciprocal basis. 




ki" 



FIG. 4: (color online). The diagram of x~^(^i5^2,fe = 0.439) 
at Tn^. It is noted from above picture the minimum point is at 
(/ci, k2) = (0.155, 0.089), which is close to the experiment value 
(0.181,0). 



effect to make ki and /c2 inequivalent. So Jnnc is of great 
importance to the magnetic groundstate, where ki is not equal 
to fe. 

Although the two layers of Mnl2 have a large separation 
about 3.5 A, it cannot be seen as a quasi-2-D triangle lattice 
of Mn atoms due to the important coupling Jnnc- 

In the mean-field approximation^ the transition temper- 
ature T/v and Curie-Weiss temperature Ocw is related to spin 
exchange parameters as 




" ki 



FIG. 5: (color online). The diagram of x~^(^i?^2 = 
0,A:3). It is noted from above picture the minimum point is at 
(/ci, ks) = (0.206, 0.444), which is close to the experiment value 
(0.181,0.439). 



TN = -^^^^J{q)min. (9) 

3kB ^ 

where the sum is over all the nearest neighbors of a given spin 
site, Zi is the number of the nearest neighbors connected by 
exchange coupling J^, and S is the spin quantum number on 
each site(in present case S = 5/2). Therefore, the T/v and 
Ocw are estimated to be 8.7 K and -10.7 K, respectively, using 
exchange parameters calculated from U=4 eV, which is com- 
pared with the experimental T/v 3.95 K and 6cw -8 K"^^. The 
ratio of the Curie- Weiss and the magnetic ordering tempera- 
ture has been proposed as a quantitative measure of frustra- 
tion. In Mnl2 the ratio is a = \Ocw /Tni | ~ 2 in experiment, 
which is comparable with that of RhFe{MoO^)2{oL ~ G)"^"^ 
and is quite a low value with respect to NaTi02(a > 500)"^^. 

To end this section, we briefly give a summary. We obtain 
the exchange parameters through DFT calculation and using 
these parameters in Heisenberg exchange model to calculate 
the magnetic vector at different magnetic phase of Mnl2 . The 
intrplane couplings lead to frustration in the triangle lattice 
of Mnl2 while the interplane coupling also has a important 
contribution. The nearest and next-nearest coupling compete 
with each other, which make the spiral magnetism stable. The 
magnetic vectors obtained are in good agreement with those 
in experiment. 



III. FERROELECTRICITY OF Mnis 
A. Symmetry analysis of Mnl2 

In this section, we perform a group theoretical calculation 
for Mnl2 magnetic structure and determine the direction of the 
FE polarization. The ferroelectric polarization appears when 
temperature is below 3 .45 K with a magnetic vector k~(0. 181, 
0, 0.439)^^. The general positions of ions with space group 
P3ml are given in Table III. 

Considering the wave vector q = q^i -\- Qzk (in our second 
Cartesian coordinate x2, y2, see Fig. 6 ), it is clear that the 
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TABLE III: General positions for P3ml. Here "3" and "2" denote 
three fold and two fold rotation, respectively. rUn label the three 
mirror plane. 





3r = {y,x + y,z) 




{x + y,x,z) 


2ir = {y,x,z) 


22Y = {x + y,y,z) 


2sr = 


{x,x + y,z) 


Ir = (x, y, z) 


I3r = (2/, X + y, z) 






miY — (y, X, z) 


m2r = {x + y,y,z) 


msr = 


= (x, X + ^, z) 



only operation (other than the identity) which conserves the 
wave vector is ms (mirror plane with respect to xz plane). 
We adopt the method in Ref.^^ to analyze the polarization of 
Mnl2. Clearly, the Fourier component Sx{(i) obeys 

msS,{q) = X{ms)S,{q), (11) 

with X{ms) = —1, and that is irrep Fi. Forirrep F2, X{ms) = 
1 and we have 

msSyiq) = X{ms)Sy{q). (12) 

To fix the complex coefficients, we consider the effect of 
inversion, which lead to 

IS^iq) = S:{q). (13) 

For Fi, we consider its quadratic free energy and substitute 
Eq.(13), then Sx{q) and ^'^(q) will have the same complex 
phase. We now introduce order parameters which describe 
the magnitude and phase of these two symmetry labels (ir- 
reps). When both irreps are present, one has 

Sx{^) = o-iiqy, 

Sz{q) = o-i{q)s, 

Sy{q) = t72(q), (14) 

where + 5^ = 1 (r,s is real) and crn(±q) = cTne^*'^"'. We 
also have the transformation properties 

msai = -t7i,m3t72 = 0-2, 

lai = al,Ia2 = cr^. (15) 

Thus, when both irreps are present, we have 

Sx{^) = cri(q)rcos((f-r + (/)i), 

Sy{r) = cr2(q)cos((f-r + (/)2), 

Sz{^) = cri(q)5cos((f-r + (/)i). (16) 

Now we consider the magnetoelectric coupling in Mnl2 us- 
ing order parameter obtained above. For single order param- 
eter cannot produce ferroelectricity in our case, we consider 
both irreps, 

F^nt = i^r^P^[^i(q)o-2*(q) - al{q)cT2{q)]. (17) 

7 

Under operation ma, cri{q)cr2{q) and crl{q)cr2{q) will 
change sign. Since that Fint is invariant under ma, it requires 
to be odd under ma, so P has to be along y direction, 
which is found in experiment^^ . 



B. Calculating the polarization using DFT 

The electronic structure of Mnl2 calculated for FM state 
with Ueff = 4 eV is presented in Fig.7. It is clear that FM 
state is insulating with an indirect band gap. Mn 3d states 
are mainly located in the lower energy region from -5.0 eV 
to -4.0 eV in the spin-up channel, while it is almost empty 
for the spin-down channel. Thus, the Mn^+ ions in Mnl2 are 
in its high- spin states. The narrow band indicates that the 3d 
electrons of Mn is highly localized. The top of the valence 
band is primarily attributed to I 5p states, hybridized weakly 
with Mn 3d, while the bottom of conduction band is mainly 
attributed to Mn 3d down spin states. The band dispersion is 
strong in ab plane but weak in c direction near fermi energy, 
as expected for layer structure of Mnl2 . 

As it has been shown that the propagation vector of Mnl2 
is q = (0.181, 0, 0.439), we perform GGA+U with the q in 
absence of SOC, in which just one primitive cell is used due to 
the generalized Bloch theroem^^. Then we calculate the elec- 
tric polarization using the Berry phase method"^^. However, 
negligible polarization is found. The above observation leads 
us to consider SOC effect on the electric polarization in the 
spiral state of Mnl2. We carry out GGA+U+SOC (Ueff = 4 
eV) calculation for the q = (0.25, 0, 0.5)"^^ spiral states with 
spin in (307) plane^"^. The polarization is 107.3 jiC /m? along 
a2 direction, which is consistent with polarization P ± q-^ 
(inplane q^^)^^. The experimental polarization is 84 jiC /w? 
at 2 K and the interpolated value at K is about 128 jiC /w?. 
The calculated polarization is a little smaller than the interpo- 
lated value, which is due to the approximation of the magnetic 
vector q. As depicted in Ref.^^, the magnetic vector is parallel 
to [100] (in our coordinate is along OD, see Fig. 6), when high 
magnetic field along [100] is applied. In this case, we perform 
calculation with q = (|, |, 0) and the ferroelectric polariza- 
tion is about 58 jiC /m? along OD, which is very close to 57 
jiC /m? at Oh = 30° in high magnetic field^^ and consistent 
with experiment P || q^^ but not with the KNB model 

Obviously, the polarization of Mnl2 is induced by SOC, 
which is also consistent with the calculation of Xiang et. al^^ 
. Let us review the spin orbit coupling. The spin-orbit part of 
Hamiltonian is given by 

Hso = XL-S = XLJ, + ^[^+^- + ^-^+]- (18) 

where A is the spin-orbit coupling constant. When the electron 
shell contains less half-filled unpaired electrons, A is positive. 
While A is negative, if there are more than half-filled electrons 
in the shell. If the shell is half filled in a ion, the spin orbit cou- 
pling vanishes because that orbit moment is zero. Considering 
Mn^+ in Mnl2, it contains five 3d electrons. Thus we expect 
spin orbit coupling will vanish in Mn. From our calculation, 
the SOC coupling on Mn is not zero although the strength is 
rather weak. However, the SOC is expected to be strong on 
/ 5p orbitals as the A increases with the nuclear charge of the 
atom and decreases with increasing quantum numbers (angu- 
lar quantum number and principle quantum number). 

To examine how the polarization arises from the spiral mag- 
netic state with SOC, we analyze the electron distribution of 
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FIG. 7: (color online). Band structure and density of states(DOS) for Mnl2 calculated with GGA+U (U=4 eV). Left panel is the band structure, 
the blue dash lines and red solid lines denote spin-up band and spin-down band, respectively. F is the center of the Brillouin zone , M (0,0.5,0), 
K( 1 , 1 ,0), A (0,0,0.5) and L (0,0.5,0.5) in reciprocal lattice. Right panel is the DOS, the positive and negative value of DOS denote the majority 
spin states and the minority spin states, respectively. 



FIG. 8: (color online). Perspective view of an isosurface calculated 
for the electron density difference between the GGA-i-U-i-SOC and 
GGA-hU results for the case of q=(0.25, 0, 0.5). The red and blue 
surface represent 5.5x 10~^e/A^ and -5.5x 10~^e/A^, respectively. 



Mnl2 by showing the difference between the electron den- 
sity of SOC calculation and that of without SOC calculation. 
For the case of k = (0.25, 0, 0.5), Fig. 8 shows that the main 
asymmetric charge distribution is around each I ion, which 
makes a primary contribution to the ferroelectric polariza- 
tion. The graphical software XCrysDen^^ was used to plot 
the electron density difference. For the case of q = (|, |, 0), 
the electron density difference is is similar to the case of 
q = (0.25,0,0.5). Furthermore, we find a linear relation- 
ship between the magnitude of electric polarization and the 
strength of SOC. 



To examine the effect of ion displacement in the spin- spiral 
state on the FE polarization, we optimize the atoms positions 
of Mnl2 in the above two cases by performing GGA+U+SOC 
calculation until the atomic forces become less than 0.02 
eV jh and then calculate the electric polarization using the 
relaxed structures. In the case of g = (0.25,0,0.5), it is 
found that Mn^+ ions move along a2 direction, which leads 
to slightly enhanced FE polarization about 170 jiC /m?. For 
q = (|, ^, 0), it is found that the sum of Mn displacements 
is along OD direction while that of I ions is along the oppo- 
site direction. The in-plane electric polarization of the relaxed 
structure is about 170 jiC /m?, but the direction of the polar- 
ization is switched. The out-of-plane component in this case 
is about one third of the in-plane component. 



IV. SUMMARY AND CONCLUSIONS 

In this paper we have presented a comprehensive investi- 
gation of the incommensurate magnetic structure and ferro- 
electric polarization of the new multiferroic material Mnl2. 
Six exchange parameters among local moments on Mn sites 
are obtained by mapping the mean-filed Heisenberg model 
Hamiltonian onto the total energy differences of eight dif- 
ferent magnetic ordering states from DFT calculations. We 
show the special inter-plane coupling Jnnc is fairly strong 
because of its linear exchange path, as suggested before by 
Wollan et aP^. Although there is a large separation between 
the Mn-layers, the lattice of Mn thus cannot be seen as a quasi 
two dimensional lattice. Moreover, the inter-plane coupling 
Jnnc is critical to generating spiral spin structure by break- 
ing the equivalence of the spin-density wavevectors along two 
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different directions, ki and k2 in the magnetic ground state. 
Our calculation also indicates Hubbard Ueff strongly affects 
the magnetic exchange couplings. It is found that Ueff ~ 4 
eV can give results quantitatively consistent with experimen- 
tal values. For example, the Currier- Weiss temperature is es- 
timated as -10.7 K when Ueff = 4 eV, which is very close to 
-8 K in the experiment measurement"^^. 

We also use both symmetry analysis and DFT calculations 
to investigate the polarization of this material. Our study re- 
veals that SOC is essential for its ferroelectric polarization. 
Both the direction and magnitude of the polarization obtained 
from DFT calculation are in good agreement with the experi- 
mental data. Charge density difference analysis shows that the 
primary asymmetric charge distribution is around I ions, due 
to their strong SOC effect. 

The isotropic Heisenberg model considered in this paper 



provides a good description of the magnetic ordering in Mnl2 . 
This suggests that the polarization induced by spiral magnetic 
ordering has no strong feedback effect on magnetic ordering. 
This result is consistent with the observed small value of po- 
larization. 
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